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Coordinate-transformation approaches to invisibility cloaking rely on the design of an anisotropic, 
spatially inhomogeneous "transformation medium" capable of suitably re-routing the energy flux 
around the region to conceal without causing any scattering in the exterior region. It is well known 
that the inherently magnetic properties of such medium limit the high-frequency scaling of practical 
"metamaterial" implementations based on subwavelength inclusions (e.g., split-ring resonators). 
Thus, for the optical range, non-magnetic implementations, based on approximate reductions of the 
constitutive parameters, have been proposed. 

In this paper, we present an alternative approach to non-magnetic coordinate-transformation 
cloaking, based on the mapping from a nearly-transparent, anisotropic and spatially inhomogeneous 
virtual domain. We show that, unlike its counterparts in the literature, our approach is amenable 
to exact analytic treatment, and that its overall performance is comparable to that of a non-ideal 
(lossy, dispersive, parameter-truncated) implementation of standard (magnetic) cloaking. 



PACS numbers; 41.20.Jb, 42.25.Fx, 42.25.Gy 

! I. INTRODUCTION AND BACKGROUND 

Invisibility of objects to an interrogating (electromag- 

I netic, acoustic, elastic, quantum) wave illumination is 
a fascinating research topic of long-standing interest, 
with a wealth of intriguing theoretical and application- 
oriented implications. For instance, in electromagnetics 
(EM) engineering, "invisible" sources, scatterers and an- 
tennas have been investigated for several decades (see, 
e.g., Refs. [llH13l3H@ for a sparse sampling). How- 

" ever, during the last few years, interest in this topic has 
gained renewed momentum, under the suggestive associ- 
ation with the "cloaking" - concept, mainly motivated by 
the rapid advances in the engineering of "metamaterials" 
with precisely controllable constitutive (e.g., anisotropy, 
spatial inhomogeneity, dispersion) properties. Among 
the most prominent approaches to (passive) invisibility 

I cloaking, it is worth recalling those based on scattering 
cancellatioD^J , coordinate transformationsi^iiiii^ii^ (ex- 

^ perimentally demonstrated at microwave^ and visible^^ 
frequencies), anomalous localized resonancesi^, inverse 
design of scattering optical elementsi^, and transmission- 

, line networks^^. The reader is referred to Ref. [3 (and 
references therein) for a recent comparative review of 
these various approaches. 

In particular, the coordinate-transformation (also re- 
ferred to as "transformation-optics" or "transformation 
EM") approachi2iii'-i2.i2., directly related to our investi- 
gation, relies on the formal invariance of Maxwell's equa- 
tions under coordinate transformations, which allows 
a preliminary design in an auxiliary curved-coordinate 
space containing a "hole," and its subsequent transla- 
tion into a conventionally flat, Cartesian space filled by 
an anisotropic and spatially inhomogeneous "transforma- 
tion medium" that suitably bends the ray trajectories so 
as to re-route the energy flux around the concealment 



region. The reader is referred to Refs. [H and [13 (and 
references therein) for a recent collection and review of 
applications and theoretical aspects. Among the most 
interesting new twists and extensions of the basic cloak- 
ing idea above, it is worth mentioning the concepts of 
EM "wormhole"^, anti-cloak22i^, carpet cloakSiiS^iS^iSl, 
cloak at a distance^^, and open cloakSS, as well as the ex- 
tensions to acoustio^i^, elastio^^, and quanturo^^ waves. 

In view of the considerable complexity of the aris- 
ing transformation media, EM modeling of such struc- 
tures typically relies on heavily numerical (finite-element) 
simulations2i. Nevertheless, at least for canonical (e.g., 
cylindrical^§i^ and spherical) geometries, analytic full- 
wave approaches, based on suitable mappings of stan- 
dard Fourier-Bessel or Mie expansions, have been devel- 
oped. Remarkably, via these approaches, it was pos- 
sible to prove analytically that in the ideal case (im- 
plying a lossless, anisotropic, spatially inhomogeneous 
transformation medium, with extreme values of the con- 
stitute parameters ranging from zero to infinity) the 
cloaking would be perfect, i.e., without any transmission 
into the concealment region and any external scattering, 
and may be, in principle, attained at any frequency - 
not necessarily in the asymptotic high-frequency regime 
that the intuitive ray-bending picture would suggest. 
However, the inherent limitations arising from the un- 
avoidable losses^^'^^, dispersion^^'^^, perturbations'^^ and 
simplifications/reductions*'' in the constitutive parame- 
ters were also pointed out. 

The first experimental verification of coordinate- 
transformation cloaking was achieved at microwave (X- 
band) frequencies^, where the involved scales allowed 
the metamaterial fabrication via low-loss metallic split- 
ring-resonator inclusions. High-frequency scaling of 
this technological solution seems to be within reach 
for the low-THz region^i, but is limited by satura- 
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tion effects^. Thus, for the visible range, an alter- 
native route has been followed, based on the use of 
non-magnetic materials^'^'^, which has led to the ex- 
perimental demonstration of an optical cloak^^. Non- 
magnetic approaches to coordinate-transformation cloak- 
ing are based on approximate reductions of the consti- 
tutive parameters that preserve the ray trajectories in- 
side the cloak shell, at the expense of destroying the per- 
fect impedance matching with the background medium, 
which can only be partially restored by suitable tweaking 
of the extra parameters available in quadratic4^i^i^ or 
higher-orde r "^^1^^ coordinate transformations. The above 
parameter reductions also prevent application of the ex- 
act analytic approaches in Refs. 35 36 Sli, and thus their 
analytic modeling is limited to asymptotic (semiclassical) 
approximations'*^ . 

In this paper, we propose a different approach to non- 
magnetic coordinate-transformation cloaking which, ac- 
knowledging the inherent limitations of realistic meta- 
material implementations, instead of applying the co- 
ordinate transformation to a perfectly transparent vir- 
tual domain (as in standard transformation optics), con- 
siders a weakly non-transparent anisotropic and spa- 
tially inhomogeneous domain. With reference to a two- 
dimensional (cylindrical) scenario, we show that, via a 
judicious choice of the constitutive parameters and the 
coordinate transformation, it is possible to achieve a 
non-magnetic transformation-medium without resorting 
to approximate reductions, thereby maintaining the ap- 
plicability of the computationally-effective and insight- 
providing exact analytic modeling (via generalization of 
the results in Refs. [35land[36l). Moreover, via suitable pa- 
rameter optimization at a given frequency, the inherently 
non-zero scattering can be minimized so that, when the 
unavoidable non-idealities (losses, dispersion, parameter 
truncations) of metamaterial implementations are taken 
into account, the overall performance becomes compara- 
ble to that of a standard (magnetic) cloak within broad 
parametric ranges. 

Accordingly, the rest of the paper is laid out as follows. 
In Sec. ini we outline the problem formulation and the 
proposed strategy. In Sec. Illli we derive the analytic 
solutions. In Sec. llVi we present some representative 
results from our parametric studies, and compare them 
with those achievable via standard (magnetic) cloaking. 
Finally, in Sec. |Vj we provide some concluding remarks 
and hints for future research. 
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FIG. 1: (Color online) Virtual-domain-medium constitutive 
parameters ej. = e'^ (black solid) and fi'^ (red dashed) in 
optimized for near-transparency at R2 — 3Ao (7 — 3.41 • 10~^, 
p = 5.41 ■ 10"*, a = 0.361, see Sec. HVA] . 



homogeneous medium immersed in vacuum. We choose 
for the relative permittivity and permeability tensors (e' 
and /i', respectively) a rather general parametric form 

featuring several degrees of freedom and yet amenable 
to analytic solution of the relevant Helmholtz equation; 
for the transverse-magnetic (TM) polarization (magnetic 
field parallel to the cylinder) of interest here, their rel- 
evant components (in the associated cylindrical coordi- 
nates r', 0, z) are given by 



4(r')=4(/)=( — 



exp 



l^'Ar)=P(r'){ — ] exp 



a 1 

i?2 



where 7 > 0, a > 0, and 

P{r') = l-p + p^, 0<p<l. 



R2 



(la) 
(lb) 

(2) 



Note that the constitutive parameters in ([T|) are always 
positive, and locally matched with vacuum at the inter- 
face r' = i?2, whereas they may exhibit singular behav- 
iors at = 0, 



II. PROBLEM FORMULATION AND 
PROPOSED STRATEGY 

A. Virtual Space: Nearly- Transparent Domain 

In the two-dimensional (2-D) scenario of interest, we 
start considering an auxiliary virtual space {x' , y' , z) fea- 
turing a circular cylindrical domain (infinite along z) 
of radius R2^ made of an anisotropic and radially in- 



lini e' 

r'-^O ' 



Xr') 



CXD, 



lim /i;(r') 



00, < 7 < 1, 
exp(-a)p, 7 = 1, 
0, 7>1. 

(3) 

We highlight that, unlike the typical transformation- 
optics framework, our virtual domain is not perfectly 
transparent. Nevertheless, the constitutive relationships 
in ID) contain three adjustable parameters (a,7,p) that 
can be optimized (see Sec. IIV Al below) for achieving a 
nearly-transparent response at a given frequency. Figure 
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(a) (b) (c) 




x-IR^ rlR^ xIR^ 



FIG. 2: (a) Ray trajectories (thin grey curves, traced along the positive x' and y' directions) in the virtual domain, with 
parameters as in Fig. [T] (b) Coordinate transformation r = g{r') vs. r' — f{r) [cf. (f9)l- (|lip ] for Ri = R2/3. (c) Ray 
trajectories in the transformed domain. The black thick circles delimit the r = Ri and r — R2 cylindrical domains. 



[T] shows an example of constitutive parameters optimized 
for an electrical size R2 = 3Ao (with Aq denoting the vac- 
uum wavelength). The log- log scale utilized highlights 
the algebraic singular behavior of the permeability for 
r' 0, but is still not able to capture the extremely 
slow (in view of the chosen parameters) divergence of 
the permittivities. For the same configuration, Fig. [2f a) 
illustrates the ray tracing, obtained as in Ref. [ij. As 
a first indication of the nearly-transparent behavior of 
this optimized parametric configuration, one can observe 
the practically straight ray trajectories, with only a slight 
bending for those passing nearby the cylinder center. 



B. Coordinate Transformation 

Following the transformation-optics 

approachiSiiiii^ii^, in order to achieve invisibility 
cloaking, we construct a cylindrical coordinate transfor- 
mation, from the virtual space {x',y',z) to the actual 
physical space (x, y, z) of interest. 



r^g{r'), r'<i?2 



(4) 



which compresses the cylindrical nearly-transparent re- 
gion r' < i?2 into a concentric annulus Ri < r < i?2, i-e., 
satisfies the boundary conditions 



3(0) <?(i?2)=i?2. 



(5) 



As previously mentioned, the non-flat metric underly- 
ing the transformation in ([4]) and ([5]) creates a "hole" 
of radius i.e., a region of space effectively impen- 
etrable by the EM fields where an object can be con- 
cealed. Invoking the formal invariance of Maxwell's equa- 
tions under coordinate transformations, the above behav- 
ior can be equivalently obtained in a globally flat space 
by filling up the transformed region Ri < r < R2 with 
an anisotropic, spatially inhomogeneous "transformation 



medium," whose relative permittivity and permeability 
tensors (e and /i, respectively) are given b yiOiHii^iis 



■l^ [det (l)] 



(6) 



where J = d{x, y, z)/d{x' , y' , z) is the Jacobian matrix of 
the transformation in ([¥]), the superfix indicates trans- 
position, and the symbol dct(-) denotes the determinant. 
For the TM polarization of interest here, the relevant 
(cylindrical) components can be directly found from ([1]) 
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.(r) = 4(r')<?(r')- 



g{r') r' ' 
l^'Ar') r' ^ 
g{r') r ' 



(7a) 
(7b) 
(7c) 



where the overdot ' denotes differentiation with respect 
to the argument, and r' = /(r) [with /(r) = g~^{r') 
denoting the inverse mapping]. In standard approaches 
to non-magnetic cloaking^^iiiii^d&iii, based on perfectly- 
transparent virtual domains (i.e., e'^. = e'^ = [i'^ = 1), 
reliance is made on the dependence of the ray trajecto- 
ries on the products ErMz ^-nd e^/Zz in order to find a set 
of reduced parameters featuring /z^ = 1. As anticipated, 
this destroys the perfect transparency of the cloak, and 
prevents application of the exact analytic framework in 
Refs. [3^ and IsgI . In our approach, instead, we capi- 
talize on the extra degrees of freedom available in order 
to enforce /^^(r) = 1 in ([7c)l . which yields the following 
first-order differential equation 



air') 



9{r') 



0. 



(8) 



A similar approach was used in Ref. '49' in order to obtain 
a spatially homogeneous axial permeability (or permit- 
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tivity, for the transverse-electric polarization). However, 
relying on a vacuum (i.e., e' = /i' = 1) virtual domain, 
such approach did not provide enough degrees of freedom. 
As a consequence, once the cloaking boundary conditions 
in ([5]) were enforced, the value of the transformed axial 
permeability was inevitably dictated by the shape fac- 
tor R1/R2, and (most important) always larger than the 
vacuum value (i.e., /iz > 1, for the case of interest here). 
In our approach, this limitation is overcome via the ex- 
tra degrees of freedom endowed by the virtual-domain 
permeability fi'z{r') [cf. Eqs. pb)) and In particu- 
lar, thanks to the functional form in (|lb[) . the differential 
equation in ([8]) admits a closed-form analytic solution as 



where 



h { — ;a,7,p 



1-P 
7 + 2 

P 



2exp(— a) 

7+2 



(9) 



M 7 + 2, 7 + 3, 



ar 
I 



r \ / ar 

7 + 1/V-R2/ V ^2 



,(10) 



with M{-, •, •) denoting the confluent hypergeometric 
function [Eq. (13.1.2) in Ref. 50]. Note that, in light 
of the assumed constraints (a > 0, 7 > 0, < p < 1), 
the function h{r' / R2; a,j,p) in pUj) is always positive, 
which ensures that the mapping g{r') in ^ is always 
real. Moreover, in © and PH)) . the arising integration 
constant has been exploited to enforce the first boundary 
condition in ([5]); the second boundary condition, instead, 
yields 



(11) 



which can be satisfied by properly tweaking the param- 
eters a, 7, and p. In particular, in view of the linear 
dependence involved, Eq. (fTTj) may be straightforwardly 
solved with respect to p. However, this seemingly sim- 
plest approach is not necessarily the most effective in a 
broader perspective of achieving a nearly-transparent re- 
sponse. In fact, our parametric studies (see Sec. IIV Al 
below) indicate that it is generally more convenient to 
satisfy the constraint in by fixing a (via numerical 
solution), and exploit the parameters 7 and p for mini- 
mizing the scattering response. 

It should be noted that, since fJ,'^{r') in (|lb[) and g{r') in 
^ are always positive, the constraint in ([5]) also implies 
g{r') > 0, and hence the invertibility of the coordinate 
mapping. However, the inverse mapping /(r) — g~^{r') 
cannot generally be calculated analytically. Neverthe- 
less, as can be observed from the example in Fig. [Hb), 
the mapping behavior is fairly regular, and therefore its 



numerical inversion does not pose any problem. Figure 
l^c) shows the ray trajectories obtained by transforming 
[via the mapping in Fig. [D(b)] those in the virtual space 
[cf. Fig. HUa)], from which the ray bending aroimd the 
interior region (r < i?i) to conceal (typical of coordinate- 
transformation-based cloaking) is fairly evident. 



C. Real Space: Non-Magnetic Transformation 
Medium 

By substituting ^ and r' = f{r) in ([7]), and taking 
into account ([5]), we readily obtain the explicit expres- 
sions of the constitutive parameters of the desired non- 
magnetic transformation medium in the real space. 



e^(r)= 



e0(r): 



Rl [fir)]'-' 



rf{r) 
Rlrf{r) 



exp 



'f{r) 



exp < —a 



R2 

R2 



(12a) 

(12b) 
(12c) 



Note that, since /(r) and f{r) are always positive in the 
transformed region Ri < r < R2, both permittivity com- 
ponents in p2ap and (|12b|l are likewise positive, thereby 
yielding a double positive medium. Moreover, like the 
virtual-space medium in ([T]), our non-magnetic transfor- 
mation medium in (|12p is locally matched with vacuum 
at r — R2. In view of the singular behavior exhibited 
by the virtual-domain medium at r' = [cf. ([3])], it is 
interesting to investigate the behavior of the transforma- 
tion medium at the image point r = Ri. Recalling that 
M(a,6,0) = 1 [cf. Eq. (13.1.2) in Ref. we obtain 
from pO| . in the limit r' 0, 

a, 7,p) - 2 exp(-a) (^4^) 



where 



0, p = 0, 

1, P^O. 



(13) 



(14) 



By substituting (I13p into ([9]), and approximating the 
square root via its first-order Taylor expansion, we then 
obtain 



g{r') ~ i?iWl 



2 cxp{-a)xRr^ (r')''-"+2 



(7 - s + 2)Rl 



Ri 



cxp(-a)xi?r''(^')''^'+^ 



(15a) 



(7 - ,s + 2)i?i 
and hence, via differentiation, 

g{r') ~ exp(-a)xi?r'^r^(^')""^+'- (15b) 
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FIG. 3: (Color online) Non-magnetic transformation-medium 
constitutive parameters Sr (black solid), (blue dotted) and 
(red dashed) in (|12p . obtained from the nearly-transparent 
virtual domain in Fig. [T]via the mapping in Fig. [2jb). 



Recalling that f(r) = l/g{r') and r' — f{r), and substi- 
tuting into (I12ap and (jl2bp . we then obtain 



Srir) 



exp(2a)i?fi?2 

X 



(16a) 



2 n27-s 



[f{r)r-^''-\ (16b) 



from which, recalling that < .s < 1 and 7 > 0, it finally 
follows 



er(i?i) =0, 



lim £^{r) 

7 — >Ri 



00. 



(17) 



Thus, as in standard (magnetic) cloaking, at the inner 
interface r = Ri, the permittivities exhibit either zero 
or infinite values. Figure [3] shows the transformation- 
medium constitutive parameters obtained from the (op- 
timized) virtual-domain medium in Fig. [1] via the coor- 
dinate mapping in Fig. [Ifb), and corresponding to the 
ray trajectories in Fig. (He). 

It it worth pointing out that, in principle, different 
choices of the parameters and functional forms of the 
virtual-domain medium ([1]) are possible, yielding finite 
variation ranges of the resulting transformation media. 
Our choice here was only aimed at facilitating direct com- 
parison with standard (magnetic) cloaks. Nevertheless, 
the effects of unavoidable parameter truncations are in- 
vestigated in our parametric studies (see Sec. IIVBI be- 
low). 



III. ANALYTIC SOLUTIONS 

The EM response of the cloaking configuration of in- 
terest can be computed analytically by generalizing the 



Fourier-Bessel-type approach in Refs. |35| and |36|. In 
what follows, we address such generalization for the 
case of TM-polarized plane-wave illumination with unit- 
amplitude z— directed magnetic field and time-harmonic 
[exp(— icjt)] dependence, by suitable coordinate mapping 
[via the inverse of ^] of the virtual-space solution. 



A. Virtual Space 

It is expedient to expand the incident plane wave H], 
(assumed to impinge from the positive x' direction) and 
the scattered field Hi into Fourier-Bessel series, 



Hlir' 



00 

i'^Jm {kor') exp(zm0), (18a) 



m— — oo 



exp{im(l)), (18b) 



m— — oo 



where ka ~ uj/c = 2Tr/ Xq denotes the vacuum wavenum- 
bcr (with c denoting the speed of light in vacuum), Jm(-) 

and H^\-) are the m-th order Bessel and Hankel func- 
tions of first kind, respectively (cf. Sec. 9.1 in Ref. 
and Cm are unknown coefficients. The magnetic field Hi 
transmitted inside the cylindrical domain r' < R2, ruled 
by the Helmholtz equation- ^ 



d 



r' n'^{r') dr' 



d 



e'^ir') dr' 
1 
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,2 ' '^o 



iI*(r',0) = O, (19) 



can be likewise expanded into this type of series 

00 

Hl{r',4>)= [a„M\\kor') + bm^^i\kor') 



m— — 00 



X exp(im0), r' < R2. (20) 

In ()20|) . ttm and bm are unknown coefficients, while 

(1 2) 

'i'm (•) denote independent solutions of the radial 
Helmholtz equations 

1 1 de'Jr') 




e;(r') dr' 



kUir')e'4r') 



m^^'4,{r') 
r''^e'^[r') 



^m(fcor') = 0. (21) 



It can be shown (see Appendix for details) that, for 
the constitutive parameters in ([1]) , these solutions can be 
expressed in closed form as 



^iYKk,r') ^ (j^J exp 



{a + ^R2)r' 



2i?2 



M{Cm,^r,r + l,^r') 
U{Cm,iym + l,^r') 



(22) 
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where (the already defined) M(-, •, •) and [/(-, •, •) are con- 
fluent hypergeometric functions [cf. Eqs. (13.1.2) and 
(13.1.3), respectively, in Ref. [HOl, and 



(23a) 



e = i?2"\/a2 - A{l-p)klRl < arg(0 < ^, (23b) 



^{v„, + l)i?2 + a(7 + 1) - 2pfc2^i 



(23c) 



Recalling that the wavefunctions exhibit a singular 
behavior at r' = (see Appendix |B|) , the field-finiteness 
condition yields 6m = in (|20p . The remaining unknown 
expansion coefficients (a„j and Cm) can be computed by 
enforcing the continuity of the magnetic and electric field 
tangential components at the interface r — R2, viz., 



i/*(i?2» = ir:(i?2,0)+i?f(i?2,0), 



(24a) 



EAR2, 



E:^{R2,(j>) + E;iR2,^), (24b) 



where the tangential electric fields readily follow from 
the curl Maxwell equation 



i?0(r-'» = 



1 



iLueos'^{r') dr' 



(25) 



By substituting the Fourier expansions (fT5|) and (PO]) into 
([M)) [with P5|) ]. we obtain a doubly countable infinity of 
linear equations, 

a,„*(;l)(fcoi?2) = t"'JmikoR2) + c,„i7£) (fcoi?2), (26a) 
a„*W(fcoi?2) = z"i,„(fcoi?2) +c™ij£)(fcoi?2), (26b) 

which can be solved in a straightforward fashion, yielding 

2jm+l 



(27a) 
(27b) 



where Ty[i?,G] represents the Wronskian of the functions 
F and G evaluated at koR2, 

W[F,G] = F (fcoi?2) G (fcoi?2) - F {koR2)G (fcoi?2) . (28) 



B. Real Space 

We can now address the solution in the real space 
(r, (j), z), by generalizing the approach in Refs. [3^ andlSGl. 
First, we note that, since the coordinate transformation 
in ([4]) is restricted within the cloak shell Ri < r < R2, 
the fields in the (vacuum) exterior region r > R2 admit 
the same expressions as in p8)) (with r substituting r'). 
For the same reason, the field transmitted into the (vac- 
uum) concealment region r < can be expanded in 
terms of a standard Fourier-Bessel series. 



dm Jm(fcor) exp(im(/)), r < Ri, (29a) 



m— — 00 



where dm are unknown coefficients, and the field- 
finiteness condition has been enforced. Following Refs. 
[35I and |3^ the field transmitted into the cloak shell 
i?i < r < i?2 is obtained by (inverse) coordinate mapping 
[via ([4])] of the virtual-space solution ([20|. viz., 



h: 



(r,0) = ^ [a„i^(^)(/cor)+6™V^)(fcor) 
Ri<r< R2, 



m— — oo 



(29b) 



where am and bm are unknown expansion coefficients, 
and 



^(^''\kor) = ^^^'^^ [fco/(r)] 



(30) 



Once again, the four sets of unknown expansion coef- 
ficients {am, bm, Cm, dm) Can be computed by enforcing 
the continuity of the tangential fields, this time at the 
interfaces r — Ri and r — R2. While the continuity at 
the outer interface [r = R2) does not pose any particular 
problem, much more involved is dealing with the inner 
interface (r = -Ri), in view of the singular behavior ex- 

hibited by the wavefunctions t/iot in (|30p [cf. Eq. (13.1.3) 
in Ref. ISO, and Appendix [B] and the azimuthal permit- 
tivity £0 in (|12b[) [c f. (fT7)) ] . In order to circumvent this 
problem, as in Refs. [ssl andfSGl. we therefore follow a lim- 
iting approach, slightly shifting the inner cloak boundary 
to r = i?! + A, where A denotes a small quantity (which 
eventually we let tend to zero). Accordingly, we obtain 



rfmJm[fco(i?i + A)] = a™VL^[fco(i?i + A)]+6^ 

rfmim [fco(i?l + A)] 



am'^m\ki){Ri 



A)] 
A)] 



^L'^[fco(i?i 

/•,(2) 



niTm 
bmlp'. 



MRi 



A)], 
-A)] 



am^lnKkoR2) + bm'4'^m\koR2) 
ami^^mHkoR2) + bmi>'i\koR2) 



i"\Jm{koR2) 
z"i„(fcoi?2) 



e4Ri + A) 

CmH^mHkoR2), 
CmH^HkoR2). 



(31a) 

(31b) 

(31c) 
(31d) 
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FIG. 4: (Color online) Curves in the (7,p) plane yielding 
a? = roots of for the parameters in Fig. [T]and shape- 
factor values R1/R2 = 1/4, 1/3, 1/2 (black-solid, red-dashed, 
blue-dotted, respectively). The regions below the curves are 
associated with positive a values, and thus represent the ad- 
missible (7, p) search spaces in the scattering- width minimiza- 
tion problem. 



It can be shown (see AppendixiBjfor details) that, in the 
limit A ^ 0, Eqs. (|3Ta|) and (|3Tb|) yield 

bm ~ amfli {-^j I (32a) 

dm ^ a,„ri2 ( ^ ) ' (32b) 

where r2i^2 are irrelevant constants. Recalling that < 
■s £ Ij 7 > and Vm > 0, it is readily realized that the 
expansion coefficients bm and dm in (j32p vanish as A — > 0. 
As a consequence, recalling (|30|) and that /(i?2) = R2, 
Eqs. picp and (|31d|) are readily recognized to become 
identical to (|26ap and (j26bp . respectively, and therefore 
the remaining unknown coefficients a„i and Cm are still 
given by (P7)) . 

To sum up, the above analytic solutions resemble those 
obtained for the standard (magnetic) cloak'^'^ in the exact 
suppression of the field transmitted into the concealment 
region (i.e., dm = 0). The expectable differences show 
up in the non-zero scattering coefficients which are 
directly inherited from the nearly-transparent medium 
^ in the virtual space. 

In what follows, we show that, via judicious optimiza- 
tion at a given frequency, the scattering response can 
be effectively reduced so that the overall performance of 
the proposed nonmagnetic cloak becomes comparable to 
that of a standard (magnetic) cloak when the unavoidable 
non-idealities (parameter truncations, losses, dispersion) 
are taken into account. 




FIG. 5: (Color online) Contour plot of the scattering width in 
(|33p as a function of 7 and p, for the parameters in Fig. [T]and 
shape factor R1/R2 = 1/3. The blue-dashed curve delimits 
the admissible search space (cf. Fig. |4j, whereas the cross 
indicates the minimum. 

IV. REPRESENTATIVE RESULTS 

We now move on to the presentation of the salient re- 
sults from an extensive series of parametric studies, start- 
ing with the virtual-domain parameter optimization, and 
continuing with the performance comparison between the 
proposed approach and standard (magnetic) cloaking. 

A. Virtual Space: Optimization for 
Near- Transparency 

As previously mentioned, the inherently non-zero scat- 
tering response of our proposed non-magnetic cloak is di- 
rectly inherited by the non-transparent virtual domain in 
([T]) . In our study, such response is compactly parameter- 
ized in terms of the total scattering cross-sectional width 
per unit length^ 

4 

Qs^TT Y] 1^™!"' (33) 

Kq 

rn— — 00 

where the scattering coefficients Cm are given by (j27b|) . 
Our approach is based on the minimization of the scatter- 
ing width in (|33p . at a given frequency, by acting on the 
free parameters (a,7,p) available in the virtual- domain- 
medium constitutive parameters ([T]). Clearly, different 
observables (e.g., more directly tied to the near-field or 
angular distributions) may be considered, giving rise to 
different optimized configurations. It is also worth recall- 
ing that the above parameters are actually constrained 
via the cloak condition in (jlip . In our numerical studies, 
we found it most effective to satisfy this constraint by 
fixing a via numerical solution of pip , and exploit the 
parameters 7 and p for minimizing Qg. Figure [3] shows. 
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FIG. 6: (Color online) Magnetic field (real part) map in the 
virtual space, for unit-amplitude TM-polarized plane-wave il- 
lumination (impinging along the positive a;— direction), and 
with ideal (lossless, non-truncated) parameters as in Fig. [T] 



FIG. 7: (Color online) As in Fig. |6l but for the real domain, 
with ideal (lossless, non-truncated) parameters as in Fig. [S] 



for representative values of the shape factor R1/R2, the 
(7,73) parametric ranges for which (|lip admits a solu- 
tion a > (consistent with the model assumptions), 
which constitute the search spaces in our minimization 
problem^-''. In view of the reduced number of parame- 
ters involved and the computationally-effective analytic 
modeling, such minimization can be readily pursued via 
exhaustive parameter scanning. For a fixed frequency 
and shape factor. Fig. [5] shows the scattering width in 
P3p as a function of the two free parameters 7 and p, 
within the above defined admissible ranges, from which 
a rather broad minimum is identified. The corresponding 
parameter configuration yields the constitutive relation- 
ships shown in Fig. [T] and the associated (practically 
straight) ray trajectories in Fig. [2ja). For the same ideal 
(lossless, non-truncated) parameter configuration. Fig. [6] 
shows the real part of the magnetic field map (for plane- 
wave excitation) computed via the expansions in (jlSp and 
([20|) [with ((27|) ]. from which the near-transparency is ev- 
ident. 



B. Real Space: Comparison with Standard 
(Magnetic) Cloak 

Referring to the non-magnetic ideal (lossless, non- 
truncated) constitutive parameters in Fig. [3] and as- 
sociated ray trajectories in Fig. ^c) [derived from the 
above optimized virtual-domain via the mapping in Fig. 
[DJb)], Fig. [7] shows the corresponding field map, com- 
puted via the expansions in (|18p and (p9)l [with (|27p and 
([32]) ]. which highlights the cloaking effect accompanied 
by a very weak scattering (identical to that in Fig. [6|). 

Although the present prototype study was not focused 
on practical applications, we did explore the effects of the 
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FIG. 8: (Color online) Parameters as in Fig. [S] but trun- 
cated at r = iii + A. Total scattering width (normal- 
ized with respect to the vacuum wavelength), as a function 
of the truncation parameter A/Ri, for loss-tangent values 
tanS = 0, 10"^10"^ (full black squares, red circles, and 
blue triangles, respectively). As a reference, the correspond- 
ing values attainable from a standard (magnetic) cloak [cf. 
(|34p ] with comparable parameter truncations are displayed 
as empty markers. 



unavoidable non-idealities, namely, the parameter trun- 
cations and material losses. As for the standard (mag- 
netic) cloak, suitable truncation of the permittivities in 
(fT2)l at the inner interface r = i?i is necessary in view of 
their singular behavior [cf. (|17p]. As in Refs. 'sH' and'36', 
in our parametric studies, we truncated the cloak shell 
at an interface r — Ri + A, considering the thin annu- 
lus Ri < r < i?i -f A as part of the concealment region. 
Figure [5] shows, for an empty (vacuum) concealment re- 
gion, the scattering width as a function of the truncation 
parameter A/i?i, for various values of the loss-tangent 





FIG. 9: (Color online) As in Fig. [H but with the conceal- 
ment region entirely filled by a dielectric object with relative 
permittivity Sobj = 4. In order to better visualize the cloak- 
ing eflFect, the scattering width is normalized with respect to 
reference value Qi°''''^ exhibited by the object in vacuum. 



ranging from zero to 10 ^. For comparison, we also stud- 
ied the response of a standard (magnetic) cloak^, 



(r) 



(ref), s 



r-Ri 

R2 



Ra — R^ 



r-Ri 



(34a) 
(34b) 

(34c) 



likewise truncated so as to guarantee comparable varia- 
tion ranges of the permittivities^^-. As one can observe, in 
spite of its non-magnetic character, the proposed cloak 
turns out to outperform the standard one (in terms of 
smaller scattering width) over a wide parametric range 
extending up to moderate values of the truncation pa- 
rameter and losses. 

It is worth pointing out that parameter truncations al- 
low field penetration through the cloak shell, and thus, 
unlike in the ideal case, the overall response generally de- 
pends on the possible presence of objects in the conceal- 
ment region. We therefore studied some representative 
scenarios, with the concealment region entirely filled up 
by a dielectric (with relative permittivity Sobj = 4 and 
16) or a perfect electric conductor (PEC), still amenable 
to analytic solution via straightforward generalization of 
([29]) . Figures [9UTT] show the corresponding responses, 
where, in order to better visualize the cloaking effect, the 
scattering width has been normalized with respect to the 
reference value qI"^"'^ exhibited by the object in vacuum. 
For all cases, it can be observed that the performance of 
proposed non-magnetic cloak is comparable to that of the 
standard cloak over wide parametric ranges. The sensible 
differences that may exist in certain parametric config- 
urations (e.g., the PEC case in Fig. [TT] for very small 



FIG. 

Sohj 



10: (Color online) As in Fig. |9l but for an object with 
= 16. 




FIG. 11: (Color online) As in Fig. [H but for a PEC object. 



values of the truncation parameter A/i?i and no losses) 
tend to disappear when moderate values of the truncation 
parameter and loss-tangent are considered. In particu- 
lar, for tan (5 = 10~^, the proposed non-magnetic cloak is 
still capable of reducing the scattering width of over an 
order of magnitude with respect to the uncloaked case, 
slightly outperforming a comparably-truncated standard 
cloak. As an example. Fig. [T^] compares the responses 
pertaining to a PEC cylinder free-standing in vacuum 
and cloaked with the proposed and standard approaches, 
for a truncation parameter A/i?i — 10~^ and loss tan- 
gent tan (5 = 10^^, in terms of the bistatic scattering 
width 



lim 27rr 
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(35) 
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FIG. 12: (Color online) Bistatic scattering width in (|35|l (nor- 
malized with respect to the vacuum wavelength), for a con- 
cealment domain filled by a PEC object. Object cloaked via 
the proposed (black-solid; Qs/Aq = 0.32) and standard (red- 
dashed; Qs/Xo — 0.41) approach, respectively, with trunca- 
tion parameter A/Ri = 10~^ and loss-tangent tanS = 10~^. 
Also shown, as a reference, is the response pertaining to 
the object free-standing in vacuum (blue-dotted; Qi°''^''/\o = 
3.46). 



which describes the angular distribution of the total scat- 
tering width in ([55]) . While the total scattering width val- 
ues for the proposed and standard cloaks turn out to be 
comparable (Qs/Aq — 0.32 and 0.41, respectively, with a 
reduction of about an order of magnitude as compared to 
the uncloaked case) , certain differences are evident in the 
angular distributions. In particular, as compared to the 
uncloaked case, our proposed cloak attains a rather uni- 
form reduction (of almost an order of magnitude, in the 
worst case), whereas the standard cloak seems to work 
much better for backscattering {(j> — 180°) but rather 
poorly in the forward {4> — 0) direction. 

Similar results, not shown here for brevity, were ob- 
tained for different frequencies and shape factors i?i/i?2- 

To sum up, when the unavoidable parameter- 
truncation and loss effects are taken into account, the 
overall performance of our proposed strategy turns out 
to be comparable to that of a standard (magnetic) cloak. 
Moreover, it should be mentioned that the theoretically 
frequency-independent behavior of this latter is severely 
restricted by the unavoidable dispersion effects in any 
practical metamaterial implementation, so that the in- 
herently narrow-band character of our proposed design 
(optimized to work at a given frequency) does not seem 
to constitute a critical limitation. 



V. CONCLUSIONS AND OUTLOOK 

In this paper, we have presented an alternative ap- 
proach to non-magnetic coordinate-transformation-based 



invisibility cloaking. Unlike other approaches in the 
literature, the proposed strategy does not rely on ap- 
proximate parameter reductions but rather on the de- 
sign, via parametric optimization, of a nearly-transparent 
anisotropic and spatially inhomogeneous virtual domain, 
and is amenable to exact analytic treatment. 

After derivation of the relevant analytic solutions, we 
have presented a body of representative parametric stud- 
ies. Our results indicate that, when the unavoidable non- 
idealities (parameter truncations, losses, dispersion) of 
typical metamaterial implementations are taken into ac- 
count, the overall performance attainable is comparable 
to that of a standard (magnetic) cloak. 

The idea underlying the proposed strategy is rather 
general and, besides the cloaking scenarios, it may open 
up interesting perspectives for other transformation- 
optics applications, such as hyper lensing^. In this frame- 
work, it should also be emphasized that the class of con- 
stitutive relationships in ([1]) represents only one example 
of nearly-transparent media amenable to analytic solu- 
tions, and exploration of further classes is certainly worth 
of interest and currently being pursued. In particular, 
configurations featuring a larger number of parameters 
may be utilized, together with more sophisticated op- 
timization strategies, in order to achieve broadband or 
multi-band responses, and/or to enforce constraints in 
the variation ranges of the constitutive parameters of the 
transformation medium. Of particular interest are also 
the "masking" application scenarios^'i^, where one is 
interested in changing the scattering signature of an ob- 
ject (e.g., making it appear larger, smaller, or of different 
shape), and thus the desired scattering response is in- 
herently non-zero. In such scenarios, the virtual-domain 
medium may offer extra degrees of freedom exploitable 
for the design of the desired scattering signature. 



APPENDIX A: PERTAINING TO EQ. (f22l) 

Particularizing the Helmholtz equation in (PT|) to the 
constitutive parameters in ([1]), we obtain 



dr'^ 



1+7 a \ d 
r' ^R^Jd? 



(7) 



*m(fcor') 

^'™(fcor') = 0, (Al) 



which, letting t = ^r' [with ^ defined in (|23b[) ]. becomes 
f- 



dr^ \^R2 
1-P)fcg I 



and, via the mapping 



l + 7^ 


-1 


r ) 


dT 


pklR2 


m 







*„4^)=0- (A2) 



— T 2 



exp 



2^R2 



^m(T), (A3) 
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reduces to 



d 



^^ + (1 + ^™-^) — -CmJ-^mM^O, (A4) 



where <;i^2 are irrelevant constants. Then, by recalhng 
([50)1 and that 



with Vm and defined in (|23a[) and (|23cp . respectively. 
Equation (jA4[) is readily recognized to be the Kummer 
equation [cf. Eq.(13.1.1) in Ref. |5(|, and admits two in- 
dependent solutions in terms of confluent hypcrgcomct- 
ric functions [cf. Eqs. (13.1.2), (13.1.3) in Ref. 50], from 
which the flnal solutions in (|^^ follow straightforwardly 
via inverse mapping. Note that, as a possible alternative 
route, the mapping 



koT\ -(2+1) 

— I = T V 2 y exp 



*m(T) (A5) 



would lead to the Whittaker equation [cf. Eq. (13.1.31) 
in Ref. M. 



APPENDIX B: PERTAINING TO THE 
APPROXIMATIONS IN (l32ll 



We start recalling that, for r' — > 0, the wavefunctions 
in (p2)) are regular, since > 1 [cf. (|23ap ] and 
M{C,^,Vm + 1,0) = 1 [cf. Eq. (13.1.2) in Ref. [H, 
whereas are singular, since [cf. Eq. (13.1.3) in Ref. 



4,(1,2) 



(kor') 



9{r') 



and substituting into (|3ip . we obtain 



A' 

am?7i 1 



(B5) 



(B6a) 



(B6b) 



dmJm(koRl) ~ ( — 



where 771^2 are other irrelevant constants. Equations (|B6|) 
can be readily solved with respect to hm and , yielding 



sin[^(i/^ + l)]r(C„Or(l-J^,„)' 



r-'^O, (Bl) 



where r(-) is the Gamma function [cf. Eq. (6.1.1) in Ref. 
[H. We then define 



A' = /(i?i + A), 
where, in view of (|15ap . 



s-7 



A ^ g(A') -R,^ exp(-a)xi?2 

' i?i(7-s + 2) ^ ' 

Thus, combining (gl]), ((BT|) and ((B2)l . we obtain 



(B2) 



(B3) 



i?2 



/ A' Y 

^71 -^m (fco^l )( I - "^1 -'^rn (fco-Rl ) 



^2Jm{koRl) - mJmikoRi) ( — 



-K2 



^2j,n{kaRi) - mJmikoRi) 



/ \ 7 ' 



i?2 



(B7a) 



(B7b) 



^(^'2)(fcoA') ^ ,1,2 



R2 



(B4a) 



(B4b) 



from which, neglecting the higher-order terms in 
(A'/i?2) and recalling (|B3[) . the approximations in 
follow straightforwardly. 
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